Reconstructing Parallel-Beam Projections

Input Image

The options in the Input image box, shown below, let you choose the projection dataset that you want to reconstruct and a reconstruction package.

Projections dataset… Lets you choose the projection dataset that you want to reconstruct.

Reconstruction package… Lets you choose a reconstruction package — RTK (Cone Beam) or TomoPy (Parallel Beam) — for tomographic data processing and image reconstruction. See Reconstructing Cone-Beam Projections for information about reconstructing cone-beam projections.

Reconstruction Engine

The options in the Reconstruction engine box let you choose a reconstruction algorithm and the acquisition parameters.

Options for parallel-beam reconstructions
	Description
Algorithm	Dragonfly provides two algorithms for tomographic reconstruction of parallel-beam projections — one is based on solving linear algebra equations and the other is a Fourier grid reconstruction algorithm. art… Algebraic reconstruction technique (art) is an iterative reconstruction technique that starts with an initial estimate of the image. Then a set of projection data is estimated from the initial estimate using a mathematical process called forward projection. The resulting projections are compared with the recorded projections and the differences between the two are used to update the estimated image. The iterative process is repeated until the differences between the calculated and measured data are smaller than a specified preselected value. For more information, refer to: Kak et al. Principles of computerized tomographic imaging. Volume 33. SIAM, 1988. gridrec… Fourier grid reconstruction algorithm. For more information, refer to: Dowd et al. Developments in synchrotron x-ray computed microtomography at the national synchrotron light source. In Proc. SPIE, Volume 3772, (224–236),1999. Rivers ML. Tomorecon: high-speed tomography reconstruction on workstations using multi-threading. In Proc. SPIE, Volume 8506, (85060U–85060U–13), 2012.
Filter	The quality of reconstructed images can be limited by several factors that result in poor spatial resolution, low contrast, and high noise levels. However, filtering can compensate for loss of detail in an image while reducing high-frequency noise. The filter chosen for a given image reconstruction task is mainly a compromise between the reduction of noise and the preservation of detail or contrast of the reproduced image. You should note that the Ramlak and Shepp filters are high pass filters that retain edges, while the Cosine, Hamming, and Hann filters are band pass filters. They are usually used to smooth images and remove extra edges from the image. The following filters, which are applied to data in the frequency domain, are available in the Filter drop-down menu: None… Filtering will not be applied. Butterworth… The Butterworth filter is a low-pass filter and is characterized by two parameters: the critical frequency which is the point at which the filter starts its roll off to zero and the order or power. Because of the ability of changing not only the critical frequency but also the steepness of the roll-off, the Butterworth filter can both smooth noise and preserve image resolution. Cosine… The Cosine filter is obtained by multiplying the Ramlak filter by a cosine function. The advantage of the Cosine and Hann filters is that they reduce image noise. A disadvantage is that they do not fully preserve edges in the image. Hamming… The low pass Hamming filter provides a high degree of smoothing and has only a single parameter to describe its shape — the cut-off frequency. The only difference with the Hann filter is on the amplitude at the cut-off frequency. Hann… The Hann ('Hanning') filter is a relatively simple low-pass filter that is very effective in reducing image noise. However, it does not preserve edges well. Parzen… The Parzen filter is another example of a low pass filter. Ramlak… The Ramlak filter is a high−pass filter that blocks out low frequencies, which can cause blurring. If there are areas in the image where the signal changes rapidly a high−pass filter sharpens the edges and gives better contrast. A disadvantage of Ramlak and other high−pass filters is that they let pass high frequencies and as a consequence also noise. Shepp… The Shepp filter belongs to the family of low pass filters and produces the least smoothing and has the highest resolution.
Min angle	Indicates the minimum angle.
Max angle	Indicates the maximum angle.
Mask ratio	The ratio between the FOV of the camera and the size of object, which is used to generate the mask.
Rotation options	Lets you enter the rotation manually.

Pre-Processing

The Pre-processing options let you choose to apply flat-field corrections, median filtering, and other corrections to improve image quality.

Pre-processing options

Pre-processing options for cone-beam projections
	Description
Flat Field Correction	If selected, lets you normalize projection data using the flat and dark field projections. In X-ray imaging, acquired projection images generally suffer from fixed-pattern noise, which is one of the limiting factors of image quality. In conventional flat field correction, projection images without the sample are acquired with and without the X-ray beam turned on, which are referred to as flat fields and dark fields respectively. Flat field… Lets you select the flat-field image taken with the X-ray beam turned on, but without the sample. Fixed pattern noise is removed with the assumption that the detector response did not change over the CT scan time. Dark reference… Lets you select the dark-field image taken before the X-ray beam was turned on.
Median	Lets you choose to apply median filtering, at a selected kernel size and shape, to help smooth noisy data.

Output Dimensions

Shows the matrix of the output dataset as the number of pixels and dimensions in the X, Y, and Z axes, as well as the voxel size.

Output dimensions

Output Dimensions

Preview

Lets you compute previews at the selected settings, as well as import the settings from a selected preview.

Preview options

Compute Preview… Click this button to generate a one-slice preview using the selected settings.

Import Inputs from Preview… Click this button to open the Import from Preview dialog. You can select any preview to reload the selected settings.

Reconstruction

Lets you reconstruct projections at the selected settings.

Reconstruct and Load… Click this button to reconstruct the selected dataset and then load the computed reconstruction into Dragonfly.